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Die folgenden Angaben sind den vom Anmelder eingereichten Unterlagen entnommen 

Prufungsantrag gem. § 44 PatG ist gestellt 

@ Verfahren zur Bildrekonstruktion fur die Computertomographie 

@ Die Erfindung betrifft ein Verfahren fur die Computerto- 
mographie, aufweisend die Verfah rensschritte, dass zur 
Abtastung eines Untersuchungsobjekts mit einem von ei- 
nem Fokus ausgehenden konusformigen Strahlenbundel 
und mit einem matrixartigen Detektorarray zum Detektie- 
ren des Strahlenbundels der Fokus relativ zu dem Unter- 
suchungsobjekt auf einer Fokusbahn um eine Systemach- 
se bewegt wird, wobei das Detektorarray der empfange- 
nen Strahlung entsprechende Ausgangsdaten in Form 
von Strahlen in Facherstrahlgeometrie liefert, dass die 
Ausgangsdaten in Richtung der zu der jeweiligen Fokus- 
position gehorigen Tangente an die Fokusbahn gefiltert 
werden und dass die gefilterten Ausgangsdaten anschlie- 
fcend zur Erzeugung wenigstens eines Schnittbildes einer 
eine Schichtdicke aufweisenden Schicht des Untersu- 
chungsobjekts, vorzugsweise unter Berucksichtigung je- 
weils aller ein und dasselbe Voxel treffender Strahlen, 
dreidimensional ruckprojiziert werden, wobei "das 
Schnittbild aus den Ausgangsdaten gewonnene Absorp- 
tionswerte der zu der Schicht des Untersuchungsobjekts 
gehorigen Voxel fur die Strahlung des Strahlenbundels 
reprasentiert. 
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Beschreibung 

[0001] Die Erfindung betrifft ein Verfahren zur Bildrekonstruktion fijr die Cqmputertomographie, aufweisend die Ver- 
fahrensschritte 

5 

a) zur Abtastung eines Untersuchungsobjekts mit einem von einem Fokus ausgehenden konusformigen Strahlen- 
biindel und mit einem matrixartigen Detektorarray zum Detektieren des Strahlenbundels wird der Fokus relativ zu 
dem Untersuchungsobjekt auf einer Fokusbahn um eine Systemachse bewegt, wobei das Detektorarray der empfan- 
genen Strahlung entsprechende Ausgangsdaten in Form von Strahlen in Facherstrahlgeometrie liefert, 
10 b) die Ausgangsdaten werden gefiltert, und 

c) die Ausgangsdaten werden an sen LieBend zur Erzeugung wenigstens eines Schnittbildes einer eine Schichtdicke 
aufweisenden Schicht des Untersuchungsobjekts ruckprojiziert, wobei das Schnittbild aus den Ausgangsdaten ge- 
wo.nnene Absorptions werte der zu der Schicht des Untersuchungsobjekts gehorigen Voxel fur die Strahlung des 
Strahlenbundels reprasentiert. 

15 

[0002] Es sind Verfahren dieser Art bekannt. unter dem Begriff "Filtered Backprojektion" bekannt geworden, die je- 
doch als 3D- Verfahren, d. h. im Zusammenhang mit einem matrixartigen Detektorarray, keine in der Praxis als ausrei- 
chend erachtete Bildqualitat erbringen, weil wegen des Kegelformigen Rontgenstrahlenbundels sogenannte "cone beam 
artifacts" auft.ret.en. 

20 [0003] AuBerdem ist an diesen Verfahren nachteilig, dass redundante Daten, wie sie bei Spiralabtastungen mit gerin- 
gem Tischvorschub durch mehrmaliges Bestrahlen ein und desselben Voxels entstehen, nicht genutzt werden, was zur 
Folge hat, dass die dem Untersuchungsobjekt verabreichte Strahlendosis nur unvollstandig zur Bildgebung genutzt wird. 
[0004] Weiter gibt. es Uberlegungen im Zusammenhang 2D- Verfahren zur Bildrekonstruktion so vorzugehen, dass vor- 
laufige Bilder in groBer Zahl durch "Filtered Backprojektion" aus Ausgangsdaten berechnet werden, die aus zur Bildre- 

25 konstruktion ansich nicht ausreichenden Abschnitten der Fokusbahn stammen, wobei die vorlaufigen Bilder erst in ei- 
nem zweiten Schritt zu einem endgultigen Schnittbild reformatiert. werden. Diese 2D- Verfahren sind fur Detektorarrays 
mit groBer Breite, d. h. groBer Erstreckung in Richtung der Systemachse, weniger brauchbar, weil dann eine extrem 
groBe Anzahl von vorlaufigen Bildern zu verarbeiten ist, was selbst bei groBer zur Verfugung stehender Rechenleist ung 
problem at isch ist. 

30 [0005] Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren der eingangs genannten Art anzugeben, das ein Stei- 
gerung der Bildqualitat ennoglicht. 

[0006] Nach der Erfindung wird diese Aufgabe gelost durch ein Verfahren mit den Merkmalen des Patentanspruchs 1 . 
[0007] Demnach erfolgt die Filterung in Richtung der zu der jeweiligen Fokusposition gehorigen Tangente an die Fo- 
kusbahn. Es hat sich gezeigt, dass mit dieser Filterrichtung eine besonders hohe Bildqualitat erzielt werden kann. Der 

35 Wahl dieser Filterrichtung liegt die Erkenntnis zugrunde, dass das beschriebene, eine hohe Bildqualitat ennoglichende, 
auf vorlaufigen Bildern basierenden 2D- Verfahren sozusagen in ein 3D- Verfahren ubergehen wiirde, wenn die der Be- 
rechnung vorlaufiger Bilder zugrundeliegenden Abschnitte der Fokusbahn so stark verkurzi wurden, dass sie nur noch 
eineeinzige Projektion umfassen wurden, deren Daten dann in Richtung der Tangente an die Fokusbahn gefiltert wurden, 
■ und dass erwartet werden konnte, dass ein solches 3D- Verfahren dann eine vergleichbar gute Bildqualitat wie das 2D- 

40 Verfahren ermoglichen wiirde. 

[0008] Das erfindungsgemaBe Verfahren lasst sich dann besonders leicht iuiplementieren, wenn vor dem Filtern eine 
Umrechnung der in Facherstrahlgeometrie gewonnenen Ausgangsdaten in Form von Strahlen P(a, p, q) in in Parallel- 
strahlgeometrie vorliegende Paralleldaten in Form von Strahlen P(0, (5, q) (azimutales "rebinning") bzw. P(0, p, q) (voll- 
standiges "rebinning", d. h. azimutales und radiates "rebinning") vorgenommen wird. Dabei sind bezugnehmend auf Fig. 

45 3 

a der Fokus winkel 
P der Facherwinkel 

q der der z-Koordinate entsprechende Zeilenindex des Detektorsystems, 
6 = a + p der Parallelfacherwinkel 
50 p = R F sin(p) die dem Strahlabstand von der Drehachse (Systemachse) entsprechende Parallelkoordinate, und 
Rp der Radius der Fokusbahn. 

[0009] Die Ruckprojektion der Paralleldaten erfolgt nach einer bevorzugten Ausfuhrungsforrri der Erfindung in der 
Weise, dass im Zuge der Ruckprojektion fur jedes Voxel (x, y, z) fur jedes 6 e [0, ft[ fur die Strahlen P(8 + krc, p, q) bzw. 
P(9 + kK y p, q), deren Projektion entlang der Systemachse durch (x, y) geht, die Summe 



55 




J) 




60 gebildet wird, wobei 

x, y, z die Koordinaten des jeweiligen Voxels sind, 

k eine ganze Zahl entsprechend der Anzahl der in die Rekonstruktion einbezogenen Halbumlaufe des Fokus, 
p die Parallelkoordinaten derjenigen Strahlen, deren Projektionen langs der Systemachse durch die Koordinaten (x, y) je- 
weiliges Voxels (x, y, z) verlaufen, 
65 p die Facherwinkel derjenigen Strahlen, deren Projektionen langs der Systemachse durch die Koordinaten (x, y) jewei- 
liges Voxels (x, y, z) verlaufen, und 

h eine die Schichtdicke der in dem erzeugten Schnittbild dargestellten Schicht des Untersuchungsobjekts bestimmende 
Gewichtungsfunktion ist, sowie 
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d eine Funktion ist, die gleich dem Abstand des jeweiligen Strahls von dem entsprechenden Voxel (x, y) ist bzw. von dem 
Abstand des jeweiligen Strahls von dem entsprechenden Voxel (x, y) abhangt. 
[0010] Die Schreibweise 



bringt dabei zum Ausdruck, dass die Sumrnenbildung wahlweise fur durch azimutales "rebinning" oder durch vollstan- 
diges "rebinning" gewonnene Strahlen erfolgen kann, wobei die Filterung tangential zur Fokusbahn im Falle des azimu- 
talen "rebinnings" einer Filterung in p-Richtung und im Falle des vollstandigen "rebinnings" einer Filterung in p-Rich- 
tung ist. . 
[0011] Infolge der Summierung sowohl uber k als auch q ist dann sichergestellt, dass alle durch ein und dasselbe Voxel 
verlaufende Strahlen berucksichtigt werden und die dem Untersuchungsobjekt zugefuhrte Strahlendosis somit vollstan- 
dig genutzt wird. 

[0012] Eine besonders bevorzugte Ausfuhrungsform sieht vor, dass zur Ruckprojektion der Paralleldaten die auf die 
Sumrne H der Gewichte h 
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normierte Summe 



H 



9 + lex, 
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gebildet wird. Diese Vorgehensweise ermoglicht. eine nochmals verbesserte Bildqualitat, da eine mogliche Uberbetonung 
von Voxcln, die von mehr Strahlen als andere Voxel getroffen werden, beseitigt wird und somit. entsprechende Artefakte 
vermieden werden. Der CT- Wert des jeweiligen Voxels wird durch Summation uber 6 gewonnen. 

[0013] Das erfindungsgernaBe Verfahren kann beispielsweise Anwendung finden, wenn. gemaB einer Variante der Er- 
rindung die Fokusbahn eine Kreisbahn ist (Tornogrammabtastung). GemkB einer bevorzugten Variante der Erfindung ist 
die Fokusbahn jedoch eine Spiralbahn, die dadurch zustande kommt, dass der Fokus auf . einer Kreisbahn tun die Systeru- 
achse bewegt wircl und gleichzeitig eine Relativbewegung zwischen Fokus und Untersuchungsobjekt in Richtung der 
Systemachse erfolgt. Auf Basis einer solchen Spiralabtastung konnen problemlos auch groBere Volumina des Untersu- 
chungsobjekts untersucht werden. 

[001.4] Im Falle einer Tornogrammabtastung gilt fur k normalerweise k = 1 oder k = 2, im Falle einer Spiralabtastung 
wird k unter Berueksichtigung der pro Vollumlauf erfolgenden Relativverschiebung in Richtung der Systemachse so ge- 
wahlt, dass der abzubildende Bereich des Untersuchungsobjekts vollstandig erfasst ist. 

[0015] Die Erfindung wird nachfolgend anhand eines in den beigefugten schematischen Zeichnungen dargestellten 
Ausfuhrungsbeispiels naher erlautert. Es zeigen: 

[0016] Fig. 1 in teils perspektivischer, teils blockschaltbildartiger Darstellung ein mehrere Zeilen von Detektoreleinen- 
ten aufweisendes CT-Gerat, 

[0017] Fig. 2 eineii Langsschnitt durch das Gerat. gemaB Fig. 1, und 
[0018] Fig. 3 ein das "rebinning" veranschaulichendes Schaubild. 

[0019] In den Fig. 1 und 2 ist ein zur Durchfuhrung des erfindungsgemaBen Verfahrens geeignetes GT-Gerat der 3. Ge- 
neration dargestellt. Dessen insgesamt mit 1 bezeichnete.Messanordnung weist. eine insgesamt mit 2 bezeichnete Ront- 
genstrahlenquelle mit. einer dieser vorgelagerten quellennahen Strahlenblende 3 (Fig. 2) und ein als flachenhaftes Array 
von mehreren Zeilen und Spalten von Detektorelementen*- eines von diesen ist in Fig. 1 mit 4 bezeichnet - ausgebildetes 
Detektorsystem 5 mit. einer diesem vorgelagerten detektornahen Strahlenblende 6 (Fig. 2) auf. In Fig. 1 sind der Uber- 
sichtlichkeit halber nur acht Zeilen von Detektorelementen 4 dargestellt, das Detektorsystem .5 weist jedoch, was in der 
Fig. 2 punktiert angedeutet ist, weitere Zeilen von Detektorelementen 4 auf. 

[0020] Die Rontgenstrahlenquelle 2 mit der Strahlenblende 3 einerseits und das Detektorsystem 5 mit der Strahlen- 
blende 6 andererseits sind in aus der Fig. 2 ersicht.licher Weise an einem Drehrahmen 7 einander derart gegenuberliegend 
angebracht, dass ein im Betrieb des CT-Gerats von der Rontgenstrahlenquelle 2 ausgehendes, durch die einstellbare 
Strahlenblende 3 eingeblendetes, pyramidenformiges Rontgenstrahlenbundel, dessen Randstrahlen mit 8 bezeichnet 
sind, auf das Detektorsystem 5 auftrifft. Dabei ist die Strahlenblende 6 dem mittels der Strahlenblende 3 eingestellten 
Querschnitt des Rbntgenstxahlenbundels entsprechend so eingestellt, dass nur derjenige Bereich des Detektorsystems 5 
freigegeben ist, der von dem Rontgenstrahlenbundel unmittelbar getroffen werden kann. Dies sind in dem in den Fig. 1 
und 2 veranschaulichten Betriebsmodus acht Zeilen von Detektorelementen 4, die im Folgenden als aktive Zeilen be- 
zeichnet werden. Die weiteren punktiert angedeuteten Zeilen sind von der Strahlenblende 6 abgedeckt und daher nicht 
aktiv. " 

[0021] Jede Zeile von Detektorelementen 4 weist eine Anzahl IC von Detektorelementen auf wobei = bis (3^ der 
Kanalindex ist und jedem Detektorelement ein Facherwinkel p k zugeordnete ist. Der Facherwinkel des mittleren Detek- 
torelementes ist gleich Null; die Facherwinkel der beiden auBersten Detektorelemente sind $ x = +p max und p K = -|} max . 
[0022] Die aktiven Zeilen L q von Detektorelementen 4 sind in Fig. 2 mit L t bis L Q bezeichnet, wobei q = 1. bis Q der 
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Zeilenindex ist, der im Falle des beschriebenen Ausfuhrungsbeispiels- sozusagen der z-Koordinate entspricht. 
[0023] Das Rontgenstrahlenbiindel weist den in Fig. 1 und 2 eingetragenen Conewinkel (p auf, bei dem es sich um den 
Offnungswinkel des Rontgenstrahlenbiindels in einer die Systemachse, Z und den Fokus F enthaltenden Ebene handelt. 
.Der Offnungswinkel des Rontgenstrahlenbiindels in einer rechtwinklig zu der Systemachse Z liegenden und den Fokus F 
5 enthaltenden Ebene (FacherdfTnungswinkel) betragt 2P max und ist in Fig. 1 eingetragen. 

[0024] Der Drehrahmen 7 kann mittels einer Antriebseinrichtung 22 um eine mit Z bezeichnete Systemachse in Rota- 
tion versetzt werden. Die Systemachse Z verliiuft parallel zu der z-Achse eines in Fig. 1 dargestellten raumlichen recht- 
winkligen Koordinatensy stems. 

[0025] Die Spalten des Detektorsys terns 5 verlaufen ebenfalls in Richtung der z-Achse, wahrend die Zeilen, deren 
10 Breite b in Richtung der z-Achse gemessen wird und beispielsweise 1 mm betragt, quer zu der Systemachse Z bzw. der z- 
Achse verlaufen. 

[0026] Um ein Untersuchungsobjekt, z. B. einen Patienten, in den Strahiengang des Rontgenstrahlenbiindel bringen zu 
konnen, ist eine Lagerungsvorrichtung 9 vorgesehen, die parallel zu der Systemachse Z, also in Richtung der z-Achse 
verschiebbar ist, und zwar derart., dass eine Synchronisation zwischen der Rotationsbewegung des Drehrahmens 7 und 
15 der Translationsbewegung der Lagerungsvorrichtung in dem Sinne vorliegt, dass das Verhaltnis von Translations- zu Ro- 
tationsgeschwindigkeit konstant ist, wobei dieses Verhaltnis einstellbar ist, indem ein gewunschter Wert fur den Vor- 
schub v der Lagerungsvorrichtung pro Umdrehung Drehrahmens gewahlt wird. 

[0027] Es kann also ein Volurnen eines auf der Lagerungsvorrichtung 9 befindlichen Untersuchungsobjekts im Zuge ei- 
ner Volumenabtastung untersucht werden, wobei die Volumenabtastung in Form einer Spiralabtastung in dem Sinne vor- 

20 genommen werden kann, dass unter gleichzeitiger Rotation der Messeinheit 1 und Translation der Lagerungsvorrichtung 
9 mittels der Messeinheit pro Umlauf der Messeinheit 1 eine Vielzahl von Projektionen aus verschiedenen Projektions- 
richtungen aufgenommen wird. Bei der Spiralabtastung bewegt sich der Fokus F der Rontgenstrahlenquelle relativ zu der 
Lagerungsvorrichtung 9 auf einer in Fig. 1 mit S bezeichneten Spiralbahn. Die Spiralabtastung muss sich in a-Richtung 
Liber wenigstens Tl + 2(5 max erstrecken, um die vollstandige Rekonstruktion eines CT-Bildes pro Zeile von Detektorele- 

25 menten zu gestattet, sie kann innerhalb der technischen Grenzen des CT-Gerats aber auch beliebig langer sein. 

[0028] Ein Volurnen des Untersuchungsobjekts kann jedoch infolge des Umstandes, dass mehrere Zeilen von Detek- 
torelementen 4 vorhanden sind, auch im Zuge einer sogenannten Tomogrammabtastung untersucht werden, bei der keine 
Relativbewegung in Richtung der z-Achse zwischen Messeinheit 1 und Lagerungsvorrichtung 9 statt.findet (v = 0). Im 
Falle der Tomogrammabtastung ist also die GroBe des untersuchten Volumens durch die Zahl der aktiven Zeilen von De- 

30 tektorelementen 4 bestimmt. Wahrend einer Tomogrammabtastung bewegt sich der Fokus F auf einer kreisformigen Fo- 
kusbahn, die in einer im Folgenden als Mittelebene bezeichneten Ebene liegt. 

1 0029] Die Tomogrammabtastung kann in Form eines Teilumlaufs oder in Form eines Vollumlaufs erlblgen, wobei der 
Teilumlauf ein Teilumlaufintervall von wenigstens n + 2f$ max (ein Halbumlauf plus Facheroffnungswinkel) umfasst, das 
eine vollstandige Rekonstruktion eines CT-Bildes gestattet, wahrend ein Vollumlauf 2n umfasst. . 
35 [0030] Die wahrend der Spiral- oder Tomogrammabtastung aus den Detektorelementen jeder aktiven Zeile des Detek- 
torsystems 5 parallel ausgelesenen, den einzelnen Projektionen P(a, p, q) in Facherstrahlgeometrie entsprechenden 
Messdaten werden in einer Datenauibereitungseinheit 10 einer DjgitaL/Analog-Wandlung unterzogen, serialisiert und an 
einen Bildrechner 11 ubertragen. 

[0031] Nach einer Vorverarbeitung der Messdaten in einer Vorverarbeitungseinheit 12 des Bildrechners 11 gelangt der 
40 resultierende Datenstrom zu einer Schnittbildrekonstruktionseinheit 13, die aus den Messdaten Schnittbilder von ge- 
wunschten Schichten des Untersuchungsobjekts nach einem noch im Einzelnen zu beschreibenden erfindungsgemaBen 
Verfahren auf Basis der "Filtered Backprojec lion" rekonstruiert. 

[0032] Die CT-Bilder setzen sich aus matrixartig zusammengesetzten Pixeln (Pixel = picture clement) zusammen, wo- 
bei die Pixel der jeweiligen Bildebene zugeordnet. sind, jedem Pixel eine CT-Zahl in Hounsfield Units (?IU) zugeordnet 

45 ist und die einzelnen Pixel entsprechend einer CT-Zahl/Grauwertskala in einem ihrer jeweiligen CT-Zahl entsprechenden 
Grauwert dargestellt werden. Dabei veranschaulicht jedes Pixel ein Voxel (Voxel = volume element) der in dem CT-Bild 
veranschaulichten Schicht des Untersuchungsobjekts. Da infolge der Mehrzeiligkeit. des Detektorsystems 5 und gegebe- 
nenfalls der Spiralabtastung Messdaten bezuglich rnehrerer Schichten des Untersuchungsobjekts gewonnen werden, ste- 
hen 3D-Daten zur Verfugung, die im Rahmen der Erfindung einer 3 D-Ruckprojekti on unterzogen werden. Als Endergeb- 

50 nis stehen 3D-Bilddaten in Form einer dreidimensionalen Matrix beispielsweise mit den Achsen x, y, z zur Verfugung, 
wobei jedes Element der Matrix einem Voxel (x, y, z) entspricht und den der zugehprigen CT-zahl entsprechenden Grau- 
wert enthalt. Diejenigen Elemente der dreidimensionalen Matrix die den gleichen x-, y-, oder z- Wert aufweisen, stellen 
dann jeweils ein planares Schnittbild der dem maBgeblichen x-, y-, oder z-Wert entsprechenden Schicht des Untersu- 
chungsobjekts dar. 

55 [0033] Die von der Schnittbildrekonstruktionseinheit 13 rekonstruierten Bilder werden auf einer an den Bildrechner 11 
angeschlossenen Anzeigeeinheit 16, z. B. einem Monitor, dargestellt. 

[0034] Die Rontgenstrahlenquelle 2, beispielsweise eine Rontgenrohre, wird von einer Generatoreinheit 17 mit. den 
notwendigen Spannungen und Stromen, beispielsweise der Rohrenspannung U, versorgt. Um diese auf die jeweils not- 
wendigen Werte einstellen zu konnen, ist der Generatoreinheit 17 eine Steuereinheit 18 mit Tastatur 19 zugeordnet, die 

60 die notwendigen Einstellungen gestattet. 

[0035] Auch die sonstige Bedienung und Steuerung des CT-Gerates erfolgt mittels der Steuereinheit 18 und der Tasta- 
tur 19, was dadurch veranschaulicht. ist, dass die Steuereinheit 18 mit dem Bildrechner 11 verbunden ist. 
[0036] Unter anderem kann die Anzahl Q der aktiven Zeilen von Detektorelementen 4 und damit die Position der 
Strahlenblenden 3 und 6 eingestellt werden, wozu die Steuereinheit 18 mit den Strahlenblenden 3 und 6 zugeordneten 

65 Versteileinheiten 20 und 21 verbunden ist. Weiter kann die Rotatiopszeit X eingestellt werden, die der Drehrahmen 7 fur 
eine vollstandige Umdrehung benotigt,.was dadurch veranschaulicht ist, dass die dem Drehrahmen 7 zugeordnete An- 
triebseinheit 22 mit der Steuereinheit 18 verbunden ist. 
' [0037] Obwohl es grundsatzlich moglich ist, das erfindungsgemaBe Verfahren auch in Facherstrahlgeometrie zu reali- 
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sieren, wird das beschriebene CT-Gerat vorzugsweise in einem Modus betrieben, in dem das erfindungsgemaBe Verfah- 
ren in ParaUelstrahlgeometrie realisiert ist. 

[0038] Demnach werden die bei der Abtastung des fur die jeweilige Untersuchung relevanten Korperbereichs des Pa- 
tienten 8 durch Spiral- oder Tomogrammabtastung in Facherstrahlgeometrie. gewonnenen Daten zunachst in an sich be- 
kannter Weise durch ein im Allgemeinen als "rebinning" bezeichnetes Verfahren in Daten in ParaUelstrahlgeometrie um- 5 
gewandelt. Diese Umwandlung beruht auf einer Umsortierung der in Facherstrahlgeometrie gewonnenen Daten derart, 
dass aus unterschiedlichen in Facherstrahlgeometrie aufgenommenen Projektionen Strahlen entnommen und zu einer 
Projektion in ParaUelstrahlgeometrie zusammengefugt werden. In ParaUelstrahlgeometrie reichen Daten aus einem In- 
tervall der der Lange TE aus, um ein vollstandiges Bild rekonstruieren zu konnem Um diese Daten gewinnen zu konnen, 
mussen nichtsdestotrotz zur Daten in Facherstrahlgeometrie aus einem Intervall der Lange n + 2p raax zur Verfugung ste- 10 
hen, " 
[0039] Fine Projektion in ParaUelstrahlgeometrie ist in Fig. 3 veranschaulicht. Demnach nehmen alle n Parallelstrah- 
len KP X bis RP^ dieser Projektion gegenuber der x-Achse des in Fig. 3 dargestellten und mit dem gemaB Fig. 1 uberein-*" 
stimmenden Koordinatensystem den Parallelfacherwinkel 8 ein. 

[0040] Anhand des in Fig. 3 in ausgezogener Linie dargestellten Parallelstrahls RP X soil im Folgenden der Ubergang 15 
von Facherstrahl- auf ParaUelstrahlgeometrie erlautert werden. 

[0041] Der Parallels trahl RP! entstammt einer fur die auf der Fokusbahn S liegende Fokusposition F] in Facherstrahl- 
geometrie gewonnenen Projektion. Der zu dieser Projektion in Facherstrahlgeometrie gehorige, durch die Drehachse 6 
und damit. die z-Achse des Koordinatensystems verlaufende Zentralstrahl RF Z1 ist in Fig. 3 ebenfalls eingetragen. Der 
Fokusposition F L entspricht der Fokuswinkel a x ; dies ist der Winkel, den die x-Achse und der Zentralstrahl RF Z1 ein- 20 
schlieBen. Gegenuber dem Zentralstrahl RF Z1 weist der Strahl RPj deri Facherwinkel p auf. Somit ist. leicht erkennbar, 
dass ftir.den Parallelfacherwinkel 0 gilt.:. 
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6 = a + p" - 

[0042] Der rechtwinklig zum jeweiligen Parallelstrahl gemessene Sl.rahlabst.and. p von der Drehachse 6 bzw. der z- 
Achse ist durch 

p = R F sin(p) • 7 

30 

gegebcn. . 

|0043] Wie anhand des in Fig. 3 in verstarkter Linie dargestellten, durch die Drehachse 6 bzw. die x-Achse verlaufen- 
den Zentralstrahls RPZ deutlich wird, handelt es sich bei diesem S trahl um den Zentralstrahl einer in Fachergeometrie fiir 
die Fokusposition FZ unter dem Fokuswinkel cx 2 aufgenommene Projektion in Facherstrahlgeometrie. Da fur den Zen- 7 
tralstrahl einer in Facherstrahlgeometrie aufgenommenen Projektion p = 0 gilt, wird deutlich, dass fur den Fall von Zen- 1 35 
tralslrahlen gilt: 

Je nachdeni, ob ein azimutales oder vollstandiges "rebinning" durchgefuhrt wird, liegen die Parallelprojektionen in der 
Form , 

P(a,p,q) - 40 

oder in der Form . 
P(B, p, q) 

. 45 

vor, wobei 

a der Fokuswinkel 

P der Facherwinkel • 
q der der ZTKoordinate entsprechende Zeilenindex des Det.ekt.orsyst.erns, 

0 = a + p der Parallelfacherwinkel . 50 

p = R F sin(P) die dem Strahlabstand von der Drehachse (Systemachse) entsprechende Paralleikoordinate, und 
R F der Radius der Fokusbahn sind. 

[0044] In einer einer ersten Ausfuhrungsform des erfindungsgemaBen Verfahrens entsprechenden, mittels der Tastatur 
19 wahlbaren ersten Betriebsart arbeitet das beschriebene CT-Gerat auf Basis von durch azimutales "rebinning" gewon- 
nenen Projektionen. Die diesen Projektionen entsprechenden Daten werden im Falle des azimutalen "rebinnings" in P~ 55 
Richtung gefiltert, d. h. jeweils in Richtung der zu der Fokusposition des Zentralstrahls der jeweiligen Parallelprojektion 
gehorigen Tangerite T (siehe Fig. 3), und zwar unter Verwendung eines der in der Computertomographie ublichen Filter- 
kerrie, z. B. Shepp-Logan- oder Ramachandran-Lakshminarayanan-Kern. 

[0045] Die derart gefilterten Paralleldaten werden dann in der in der Weise riickprojiziert., dass im Zuge der Riickpro- 
jektion fur jedes Voxel (x, y, z) fiir jedes 9 e [0, n[ fur die Strahlen P(8 + kTC, P, q), deren Projektion en tlang der System- 60 
achse durch (x, y) gent, die Summe 

k q • - ■ 

gebildet wird, wobei 

x, y, z die Koordinaten des jeweiligen Voxels, sihd, 

k eine ganze Zahl entsprechend der Anzahl der in die Rekonstruktion einbezogenen Halbumlaufe des Fokus, 
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P die Facherwinkel derjenigen Strahlen sind, deren Projektionen langs der Systemachse durch die Koordinaten (x, y) je- 
weiliges Voxels (x, y, z) verlaufen, und 

h eine die Schichtdicke der in dem erzeugten Schnittbild dargestellten.Schicht.des Untersuchungsobjekts bestimmende 
Gewichtungsfunktion ist, sowie 
5 d eine Funktion ist, die gleich dem Abstand des jeweiligen Strahls von dem entsprechenden Voxel (x, y) ist bzw. von dem 
Abstand des jeweiligen Strahls von dem entsprechenden Voxel (x, y) abhangt. 

[0046] Infolge der gewahlten Filterrichtung und der infolge der Summation sowohl uber k als auch uber q werden zum 
einen cone beam artifacts vermieden und zum anderen im Interesse einer hohen Dosisnutzung aller durch ein Voxel (x, y, 
z) verlaufenden Strahlen beriicksichtigt. 
10 [0047] Der einem Voxel x, y, z zugeordnete Absorptionswert u x>y , z wird durch Summation uber 6 iiber mindestens ei- 
nen Halbumlauf gewonnen, d. h. durch Bildung von 

9 

15 

[0048] Die dem jeweiligen Absorptionswert entsprechende CT-Zahl wird in herkommlicher Weise aus dem Absorpti- 
onswert ermittelt. 

[0049] Dabei konnen mittels derTastatur 19 unterschiedliche Gewichtungsfunktionen h und unterschiedliche Funktio- 
nen d eingestellt werden. 

20 [0050] Als Gewichtungsfunktion h eignet sich beispielsweise eine Dreiecks- oder Trapezfunktion. 

[0051] Als Funktion d kann der Abstand des jeweiligen Parallelstrahls vom Voxel x, y, z oder stat.tdessen beispiels- 
weise die z(axiale)-Komponente dieses Abstandes eingestellt werden. 

[0052] In einer Abwandlung der vorstehend beschriebenen ersten Betriebsart wird zur Ruckprojektion der Parallelda- 
, ten die auf die Summe H der Gewichte h 

25 

H = ^^h{d x ^Xo + kn,p y q)) 

.kg 
nonnierte Summe 

30 

Ky, (*) = 77 S 2 h ( d , y , i 0+k *' P^)W e + kxJ,q) 

"kg 

gebildet wird. Dies ermoglicht eine nochmals verbesserte Bildqualitat, da eine inogliche Uberbeionung von Voxeln, die 
35 in mehreren Halbumlaufen "beleuchtet", also von Strahlen getroffen werden, beseitigt wird und somit. entsprechende Ar- 
tefakte vermieden werden. Diese Redundanz trilt bei Spiralabtastungen dann auf, wenn die pro Vollumlauf der Messan- 
ordnung stattfindende Relativverschiebung so gering ist (geringer Pitch), dass Voxel uiehrfach bestrahlt werden. 
[0053] Eine einer weiteren Ausfiihrungsfonn des erfindungsgemaRen Verl'ahrens entsprechenden, mittels derTastatur 
19 wahlbaren zweite Betriebsart unterscheidet sich von der ersten Betriebsart dadurch, dass das beschriebene CT-Gerat 
40 nicht auf Basis von durch azimutales, sondern von durch vollstandiges "rebinning" gcwonnenen Projektionen arbeitet. 
Die diesen Projektionen entsprechenden Daten werden im Falle des vollstandigen "rebinnings" in p-Richtung gefiltert, 
also jeweils ebenfalls in Richtung der zu der Fokusposition des Zentralstrahls der jeweiligen Parallelprojektion gehori- 
gen Tangente T (siehe Fig. 3). 

[0054] Demnach wird'fiir die derart gefilterten Paralleldaten im Zuge der Ruckprojektion die Summe 

k g 

gebildet, wobei 

50 p die ParallelkoQrdinaten derjenigen Strahlen sind, deren Projektionen langs der Systemachse durch die Koordinaten (x, 
y) jeweiliges Voxels (x, y, z) verlaufen. 

[0055] Auch im Falle der zweiten Betriebsart wird in einer Abwandlung zur Ruckprojektion der Paralleldaten eine auf 
die Summe H der Gewichte h 

kg 

normierte Summe, namlich die Summe 

60 P,yA^)-^^h(d x ^^kn,p, q ))p{e + kK,p,q) 

"kg 

gebildet. 

[0056] Im Falle der vorstehend beschriebenen ersten und zweiten Betriebsarten wird eine Funktions weise des erfin- 
65 dungsgemaBen Verfahreris vorgesehen, bei dem im Zusammenhang mit einem Voxel (x, y, z) alle Strahlen betrachtet 
werden, deren Projektion entlang der Drehachse 6 bzw. der ZrAchse durch x, y geht. Ob und in welchem MaBe diese 
Strahlen berucksichtigt werden, wird durch die Gewichtungsfunktion h und die Funktion d bestimmt. 
[0057] Das CT-Gerat kann aber auch mittels der Tastatur 19 wahibare weitere Betriebsarten aufweisen, die den zuvor 
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beschriebeneri mit dem Unterschied entsprechen, dass fur eine gegebene Fokusposition der durch das jeweilige Vbxel (x, 
y, z) verlaufende theoretische Strahl ermittelt wird und dann unter Beriicksichtigung der Gewichtungsfunktion h und der 
Funktion d in die Summenbildung im Zuge der Riickprojektion nur diejenigen Strahlen einbezogen werden, die tatsach- 
lich einen Beitrag zu der Summe liefem konnen. 

[0058] Im Falle der beschriebenen Ausfuhrungsbeispiele wird die Relativbewegung zwischen der Messeinheit 1 und 5 
Lagerungsvorrichtung 9 jeweils dadurch erzeugt, dass die Lagerungsvorrichtung 9 verschoben wird. Es besteht im Rah- 
men der Erfindung jedoch auch die Moglichkeit, die Lagerungsvorrichtung 9 ortsfest zu lassen und statt dessen die Mess- 
einheit 1 zu verschieben. AuBerdem besteht im Rahmen der Erfindung die Moglichkeit, die notwendige Relativbewe- 
gung durch Verschiebung sowohl der Messeinheit 1 als auch der Lagerungsvorrichtung 9 zu erzeugen. 

[0059] Das konusformige Rontgenstrahlenbiindel weist im Falle des beschriebenen Ausfuhrungsbei spiels einen recht- 10 
eckigen Querschnitt auf. Im Rahmen der Erfindung sind jedoch auch andere Querschnittsgeometrien moglich. 
[0060]. Im Zusammenhang mit den vorstehend beschriebenen Ausfuhrungsbeispielen finden CT-Gerate der dritten Ge- 
neration Verwendung, d. h. die Rontgenstrahlenquelle und das Detektorsystem werden wahrend der Biiderzeugung ge- 
meinsam um die Systemachse verlagert. Die Erfindung kann aber auch im Zusammenhang mit CT-Geraten der vierteh 
Generation, bei denen nur die Rontgenstrahlenquelle um die Systemachse verlagert wird und mit einem feststehenden is 
Detektorring zusammenwirkt, Verwendung finden, sofern es sich bei dem Detektorsystem um ein mehrzeiliges Array 
von Detektorelernenten handelt. 

[0061] Auch bei CT-Geraten der funften Generation, d. h. CT-Geraten, bei denen die Rontgenstrahlung nicht nur von 
einem Fqkus- sondern von niehreren Foken einer oder mehrerer um die Systemachse verlagerter Rontgenstrahlenquellen 
ausgeht, kann das erfindungsgemaBe Verfahren Verwendung finden, sofern das Detektorsystem ein- mehrzeiliges Array 20 
von Detektorelernenten aufweist. 

[0062] Die im Zusammenhang mit den vorstehend beschriebenen Ausfuhrungsbeispielen verwendeten CT-Gerate wei- 
sen ein Detektorsystem mit nach Art einer orthogonalen Matrix angeordneten Detektorelernenten auf. Die Erfindung 
kann aber auch im Zusammenhang mit CT-Geraten Verwendung finden, deren Detektorsystem in einer anderen Weise 
flachenhaft.es Array angeordnete Detektorelemente aufweist. 25 
[0063] Die vorstehend beschriebenen Ausfuhrungsbeispiele betreffen die medizinische Anwendung des erfindungsge- 
maBen Verfahren s. Die Erfindung kann jedoch auch aufierhalb der Medizin, bei spiels weise bei der Gepackprufung oder 
bei der Materialuntersuchung, . Anwendung finden. 



Patentanspruche j ( ) 

1 . Verfahren zur Bildrekonsiriiktion fur die Computertomographie, aufweisend die Verfahrensschritte: 

a) zur Abtastung eines Untersuchungsobjekts mit einem von einem Fokus ausgehenden konusformigen Strah- 
lenbundel und mit einem matrixartigen Dctektorarray zum Detektieren des Strahlenbundels wird der Fokus re- 
lativ zu dem Untersuchungsobjekt auf einer Fokusbahn urn eine Systemachse bewegt/wobei das Detektorarray 35 
der empfangenen Strahlung entsprechende Ausgangsdaten in Form von Strahlen in Fachers.trahlgeornet.rie lie- 
rert, 

b) die Ausgangsdaten werden in Richtung der zu der jeweiligen Fokusposition gehorigen Tangente an die Fo- 
kusbahn gefilt.ert, und 

c) die gefilterten Ausgangsdaten werden anschlieBend zur Erzeugung wenigstens eines Schnittbildes einer 40 
. eine Schichtdicke aufweisenden Schicht. des Untersuchungsobjekts dreidimensional ruckprojiziert, wobei das- 

Schnittbild aus den Ausgangsdaten gewonnene Absorptionswert.e der zu der Schicht des Untersuchungsob- 
jekts gehorigen Voxel fur die Strahlung des Strahlenbundels reprasentiert. 

2. Verfahren nach Anspruch 1, bei dem vor dem Filtern eine Umrechnung der in' Facherstrahlgeometrie gewonne- 

nen Ausgangsdaten in Form von Strahlen P(a, p, q) in in Parallels trahlgeome trie vorliegende Paralleldaten in Fonn 45. 
von Strahlen P(6, P, q) bzw. P(9, p, q) vorgenommen wird, wobei 
a der Fokuswinkel 
P der Facherwinkel 

q der der z-Koordinate entsprechende Zeilenindex des Detekt.orsyst.ems, 

6 = a + p der Parallelfacherwinkel 50 
p = Rpsin(p) die dem Strahlabstand von der Drehachse (Systemachse) entsprechende Parallelkoordinate, und 
R F der Radius der Fokusbahn sind. 

3. Verfahren nach Anspruch 2, bei dem die Riickprojektion der Paralleldaten in der Weise erfolgt, dass im Zuge der 
Riickprojektion fur jedes Voxel (x, y, z) fur jedes 9 e [0, n[ fur die Strahlen P(0.+ kK, P, q) bzw. P(B + kTC, p, q), de- 
ren Projektion entlang der Systemachse durch (x, y) geht, die Summe 55 

k q 

60 

gebildet wird, wobei 

x, y, z die Koordinaten des jeweiligen Voxels sind, 

k eine ganze Zahl entsprechend der Anzahl der in die Rekonstruktion einbezogenen Halbumlaufe des Fokus, 
p die Parallelkoordinaten derjenigen Strahlen sind, deren Projektioneri langs der Systemachse durch die Koordina- 
ten (x, y) jeweiliges Voxels (x, y, z) verlaufen, • 65 
P die Facherwinkel derjenigen Strahlen sind, deiren Projektionen langs der Systemachse durch die Koordinaten (x, 
y) jeweiliges Voxels (x, y, z) verlaufen, und 

h eine die Schichtdicke der in dem erzeugten Schnittbild dargestellten Schicht des Untersuchungsobjekts bestim- 
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mende Gewichturtgsfunktion ist, sowie 

d eine Funktion ist, die gleich dem Abstand des jeweiligen StrahLs von dem entsprechenden Voxel (x, y) ist bzw. von 
dem Abstand des jeweiligen Strahls von dem entsprechenden Voxel (x, y).abhangt. 

4. Verfahren nach Anspruch 3, bei welchem bei der Ruckprojektion der Paralleldaten die auf die Sumrne H der Ge- 
5 wichte h 



10 



15 




gebildet wird. 

5. Verfahren nach einem der Anspriiche 1 bis 4, bei dem die Fokusbahn eine Kreisbahn ist. 

6. Verfahren nach einem der Anspriiche 1 bis 4, bei dem die Fokusbahn eine Spiralbahn ist, die dadurch zustande 
20 kommt, dass der Fokus auf einer Kreisbahn urn die Systemachse bewegt wird und gleichzeitig eine Relativbewe- 

gung zwischen Fokus und Untersuchungsobjekt in Richtung der Systemachse erfolgt. 
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